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Abstract 

A non-Hermitian Hamiltonian has a real positive spectrum and exhibits unitary time evolution 
if the Hamiltonian possesses an unbroken VT (space-time reflection) symmetry. The proof of 
unitarity requires the construction of a linear operator called C. It is shown here that C is the 
complex extension of the intrinsic parity operator and that the C operator transforms under the 
Lorentz group as a scalar. 
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I. INTRODUCTION 



If the Hamiltonian that defines a quantum theory is Hermitian H = H>, then we can 
be sure that the energy spectrum of the Hamiltonian is real and that the time evolution 
operator U = e lHt is unitary (probability preserving). (The symbol f represents conventional 
Dirac Hermitian conjugation; that is, combined complex conjugation and matrix transpose.) 
However, a non-Hermitian Hamiltonian can also have these desired properties. For example, 
the quantum mechanical Hamiltonians 

H = p 2 + x 2 (ix) e (e>0) (1) 

have been shown to have real positive discrete spectra and to exhibit unitary time 

evolution [3J. Spectral positivity and unitarity follow from a crucial symmetry property of 
the Hamiltonian; namely, VT symmetry. (The linear operator V represents parity reflection 
and the anti-linear operator T represents time reversal.) To summarize, if a Hamiltonian 
has an unbroken VT symmetry, then the energy levels are real and the theory is unitary. 

The proof of unitarity requires the construction of a new linear operator called C, which in 
quantum mechanics is a sum over the eigenstates of the Hamiltonian |3]. The C operator is 
then used to define the inner product for the Hilbert space of state vectors: (A\B) = A cvr -B. 
With respect to this inner product, the time evolution of the theory is unitary. 

In the context of quantum mechanics, there have been a number of papers published on 
the calculation of C. A perturbative calculation of C in powers of e for the cubic Hamiltonian 

H = p 2 + x 2 + iex 3 (2) 

was performed 4| and a perturbative calculation of C for the two- and three-degree-of- 
freedom Hamiltonians 

H = p 2 + q 2 + x 2 + y 2 + iexy 2 (3) 

and 

H = p 2 + q 2 + r 2 + x 2 + y 2 + z 2 + iexyz (4) 

was also performed P|. A nonperturbative WKB calculation of C for the quantum- 
mechanical Hamiltonians in was done 0]. 

In quantum mechanics it was found that the C operator has the general form 

C = e Q V, (5) 



where Q is a function of the dynamical variables (the coordinate and momentum operators). 
The simplest way to calculate C in non-Hermitian quantum mechanics is to solve the three 
operator equations that define C: 

C 2 = l, [C,VT]=0, [C,H]=0. (6) 

Following successful calculations of the C operator in quantum mechanics, some calcu- 
lations of the C operator in field-theoretic models were done. A leading-order perturbative 
calculation of C for the self-interacting cubic theory whose Hamiltonian density is 

W(x, t) = |vr 2 (x, t) + !/iV(x, t) + f pv>(x, t)f + ^ 3 (x, t) (7) 

was performed [3, 0] and the C operator was also calculated for a "PT-symmetric version of 
quantum electrodynamics j^. For each of these field theory calculations, it was assumed 
that the form for the C operator is that given in (jSJ), where V is now the field-theoretic 
version of the parity reflection operator; that is, for a scalar field 

V(p(x,t)V = (p(-x,t), (8) 

and for a pseudoscalar field 

V<p{x,t)V = -<p{-x,t). (9) 

However, in a recent investigation of the Lee model in which we calculated the C operator 
exactly [10], we found that the correct field-theoretic form for the C operator is not C = e®V 
but rather 

C = SVi. (10) 

Here, Vi is the intrinsic parity operator; Vj has the same effect on the fields as V except 
that it does not change the sign of the spatial argument of the fields. Thus, for a scalar field 

Pjp(x,*)Pj = p(x,t), (11) 

and for a pseudoscalar field 

^(x,^ = -^ X) t). (12) 

The fundamental difference between the conventional parity operator V and the intrinsic 
parity operator Vj is seen in their Lorentz transformation properties. For a quantum field 
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theory that has parity symmetry the intrinsic parity operator Vi is a Lorentz scalar because 
Vi commutes with the generators of the homogeneous Lorentz group: 



[Pi, J^) = 0. (13) 
However, the conventional parity operator V does not commute with a Lorentz boost, 

[p, J 0i ] = -2J 0i V, (14) 
and V transforms as an infinite-dimensional reducible representation of the Lorentz group 



111 ]. Specifically, V transforms as an infinite direct sum of finite-dimensional tensors: 

(0, 1) © (0, 3) © (0, 5) © (0, 7) © ■ ■ • . (15) 

That is, V transforms as a scalar plus the spin-0 component of a two-index symmetric 
traceless tensor plus the spin-0 component of a four-index symmetric traceless tensor plus 
the spin-0 component of a six-index symmetric traceless tensor, and so on. Note that in 
(1151) we use the notation of Ref. 12 1. It was shown in Ref. 11] that the decomposition in 

nn 

(J15|) corresponds to a completeness summation over Wilson polynomials [111 Il3||. 

We believe that the correct way to represent the C operator is in (jlOj) and not in (|5jl. In 
the case of quantum mechanics there is, of course, no difference between these two repre- 
sentations because in this case V = Vi- However, in quantum field theory, where P ^ 
these two representations are different. For the case of the Lee model, (0) is the wrong 
representation for C and ()10|) is the correct representation. It is most remarkable that for 
the case of the cubic quantum field theory in ((7|) in either representation the functional Q 
is exactly the same. However, the representation of C in (|10|) is strongly preferred because, 
as we will show in this paper, it transforms as a Lorentz scalar. 

The work in this paper indicates for the first time the physical and mathematical in- 
terpretation of the C operator: The C operator is the complex analytic continuation of the 
intrinsic parity operator. To understand this remark, consider the Hamiltonian H associated 
with the Hamiltonian density 7i in (J7|), H = J c?x7Y(x, t). When e = 0, H commutes with 
Vi and Vi transforms as a Lorentz scalar. When e 7^ 0, H does not commute with Vi, 
but H does commute with C. Moreover, we will see that C transforms as a Lorentz scalar. 
Furthermore, since Q — ► as e — > 0, we see that C — > Vi in this limit. Therefore, we can 
interpret the C operator as the complex extension of the intrinsic parity operator. 



This paper is organized very simply: In Sec. |H] we examine an exactly solvable VT- 
symmetric model quantum field theory and we calculate the C operator exactly. We show 
that this C operator transforms as a Lorentz scalar. Next, in Sec. II 111 we examine the cubic 
quantum theory whose Hauultoniau density is gl ven in Q and show that the C operator, 
which was calculated to leading order in e in Ref. |7J, again transforms as a Lorentz scalar. 
We make some concluding remarks in Sec. IIVI 



II. THE C OPERATOR FOR A TOY MODEL QUANTUM FIELD THEORY 

In this section we calculate the C operator for an exactly solvable non-Hermitian VT- 
symmetric quantum field theory and show that C transforms as a Lorentz scalar. Consider 
the theory defined by the Hamiltonian density 

W(x,t) = W (x,t) + efti(x,t), (16) 

where 

WoCM) = i7r 2 (x,t) + i/iV 2 (x,t) + |[V x ^(x,t)] 2 , Wi(x,t)=i^(x,t). (17) 

We assume that the field (p(pc,t) in (JT7j) is a pseudoscalar. Thus, under intrinsic parity 
reflection <f(pc,t) transforms as in ()12|) . Also, 7r(x, t) = (p(x,t), which is the field that is 
dynamically conjugate to <p(x, t), transforms as 

P/7T(x,t)Pi = -7T(x,t). (18) 

To calculate the C operator we assume that C has the form in tjlOj). The operator Q[tt, (p], 
which is a functional of the fields <p and it, can be expressed as a series in powers of e: 

Q = eQ x + e 3 Q 3 + e 5 Q 5 + ■■■ . (19) 

Note that only odd powers appear in the expansion of Q. As shown in Ref. j^, we obtain 
Qi by solving the operator equation 

Q u Jdx.H {x,t) =2 y*dxWi(x,t). (20) 



To solve (|20|) we recall from Ref. (j| that the functional Q[7r, is odd in 7r and even in 
ip. We then substitute the elementary ansatz 



Qi 



fdxR(x)ir(x,t), (21) 



where -R(x) is a c-number. We find that -R(x) satisfies the functional equation 

y*dx^(x,0(/i 2 -V x 2 )i2(x) =-2 Jdbc<p{x,t), (22) 



whose solution is a constant: 

i?(x) = -2//i 2 . 

All higher-order contributions to Q[k, (p] in (|19j) vanish because 

[Qi, [Qi, H x }} = [Q 1 ,-2y/ / / 2 ] =o, 



(23) 



(24) 



where V is the volume of the space. Therefore, the sequence of equations in Eq. (34) of 
Ref. Q terminates after the first equation. We thus obtain the exact result 



C = exp ( 1 J da. 7r(x, t)^ Vi. 



(25) 



The C operator is clearly a rotational scalar. To prove that C transforms as a Lorentz 
scalar, we must show that C commutes with the Lorentz boost operator 



where 



J 0i = J Q 0i + 



J 0l (t) = t /dx7r(x,t)VXx,t) - dxx l H (x,t) 



(26) 



jf(t) = -yc/xx^i(x,t). 

The commutator of C = e^Vi with J°* in ()2fi)l has two terms: 
[C, J 04 ] = [e Q Vj, J° l + eJf ] 

= -2ee Q J Q l l V I + [e Q , J 0i + eJf] Vj. 

To evaluate the second term in (J28)) we use the general formula 
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(27) 



f[ir],(p(it,t) 



<$7t(x, t) 



(28) 



(29) 



from which we obtain 



7^ exp (-^/ rfx7r(x ' t} )- 



(30) 



Thus, the second term in (|28|) evaluates to 

[e« J °* + eJ«] = f |rf X7 r(x,t)V^e« - |dxx* Q/^x,*)^ + ^ 2 ^e^(x,t) 

+-V X (^ x, t • V x ^e Q + -V x ^e Q ■ V x ^(x, t) + ie^e Q 

= -«/ rf xxV(x, t )e«- i£ e«/ d x lV (x, t ) 

= -2iee Q Jdxx^x^t). (31) 

In the calculation above we assume that the integral J dx vanishes by oddness, and in the 
last step we use the identity 

Jdx.x l ip(yi,t), Jdy 7r(y, t) = z Jjdxdy x % 5(x — y) = z Jdxx t = 0. (32) 

Plugging (}3*T]) into the commutator (J2EJ), we obtain 

[C,J 0l ] = 0, (33) 

which verifies that the C operator for this toy model transforms as a Lorentz scalar. 

III. THE C OPERATOR FOR AN i(^ 3 QUANTUM FIELD THEORY 

We consider next the cubic quantum field theory described by the Hamiltonian density 
in (J2J) and we rewrite H as TL(x, t) = 7io( x , t) + e7ii(x, t), where 

Wo(x,t) = |vr 2 (x,t) + |/iV(x,t) + |[V x </?(x,t)] 2 , Wi(x,t) = ^ 3 (x,t). (34) 

For this nontrivial field theory the C operator still has the form C = e®Vi, where Q is a series 
in odd powers of e, Q = eQi + e 3 Qs + e 5 Q$ + • • • , and Vi is the intrinsic parity operator. 

We expect that the operator C is a rotational scalar because neither Q nor Vi depend on 
spatial coordinates and one can verify that C is indeed a rotational scalar by showing that 
C commutes with the generator of spatial rotations J* J . To prove that C is a Lorentz scalar, 
we must show that C also commutes with the Lorentz boost operator J° l = J® 1 + eJ° l , where 
the general formulas for J® 1 and Jx'it) are given in ([27)1 . 

Let us expand the commutator [C, J 0i ] in powers of e. To order e 2 , we have 

[C, J *] = [(l + eQ x + ±e 2 Q 2 ) Vi, J °* + eJf] + 0(e 3 ) 

= [Vi, J 0i ] + e ( [Qi, J 0i ] + Qi [P/, J°o l ] + [Pj, J , 1 } ) 

+6 2 (i [Q 2 , J«] ^ + \Q\ [Vi, J 0i ] + [Q X , + Qi [Pi, fl) + 0(e 3 ). (35) 



The leading term vanishes because Vi commutes with J® 1 : [Pi, J® 1 } = 0. Using the identity 
fp 7) JK] = -2JfV u we simplify flUSJ) to 

[C,J 0i ] = e([Q 1 ,J° i ]V I -2J° i T I ) 

+e 2 (|Qi [Qi, J 0i ] Pi + I [Qi, <5iPi " Qi^Pi - JTQiVi) + 0(e 3 ). (36) 

Evidently, if the term of order e vanishes, then the e 2 term vanishes automatically. 

From (J36j) we can see that to prove that C is a scalar we need to show that [Qi, J® 1 ] = 2 J® 1 . 
The generator J® 1 in (|27|) consists of two parts, and the first part commutes with Q\. To 
show that this is so, we argue that for any functional of tt and ip, say f[ir, ip], the first term 
of Jq 1 in (}2*T|) commutes with f\n,tp\: [f[ir, <p],t f g?X7t(x, t)V^cp(x., £)] = 0. We can verify 
this either by using time-translation invariance and setting t = or by noting that this 
commutator is explicitly an integral of a total derivative: 



f[ir,<p],t J dy*(y,t)V*<p(y,t) 
t dy li 



!) /[7r,^]V^(y,t)-m(y,t)V; " 



My>*) 



My,*) 



(37) 



where we have used the variational formulas [/ [tt, ip] , </?( x , t 



'Sffej/K*'] and 



' <S<p(x,t) 



f[ir, ip]. Integrating the second term of (|37j) by parts, we get 



/[7T,^],t / dy7r(y,t)V>(y,t) 



0. 



(38) 



Thus, since Qi is a functional of 7r and <p, only the second term of J® 1 in (|27j) contributes to 
the commutator of Qi with J® 1 : [Qi, J® 1 ] = — [Qi, J dxx l Ho(x,t)]. 

Thus, we have reduced the problem of showing that C is a scalar to establishing the 
commutator identity 

Qi, / dxa;^|7r 2 (x,t) + ±/iV(x,t) + |[V x ^(x, t)] 2 )l = 2i j rfxxV 3 (x,t). (39) 

T his i U a ti on issim „ Mi n_e t o Eq .( 66 , i „Ke,Q, 

g 1 ,y'rfx(|7r 2 (x,t) + | /U V(x,t) + |[V x ^(x,t)] 2 )| = 2* Jdx<p 3 (x,t), (40) 

apart from an integration by parts. The only difference between (}3*9"j) and (|40|) is that there 
are extra factors of in the integrands of ()39|) . 



We will now follow the same line of analysis used in Ref. Q| to solve (J3HJ). We introduce 
the identical ansatz for Qi. 



Qi = jjjdxdydz M( xyz )7T x 7r y 7r z + jjjdxdydz iV x ( yz )y? y 7i- x y2 z , (41) 

where we have suppressed the time variable t in the fields and have indicated spatial depen- 
dences with subscripts. To indicate that the unknown function M is totally symmetric in 
its three arguments, we use the notation M( xyz ) and we write iV x ( yz ) because the unknown 
function N is symmetric under the interchange of the second and third arguments. 
Performing the commutator in (}3*9*j) . we obtain two functional equations: 



dxdydz (p x (p y (p z [{x l + y i + z^fj? - V* - ^V x ] N x(yz) = -2 J dwufip 



w' 



(42) 



dltdydz (yV y 7T x V9 z + ^Vy^xTTz) iV x ( yz ) 



dx dy dz 



V? x 7r y 7r z (x> 2 - V* - a^V 2 ) M (xyz) 
+^ y n z (y> 2 - V* - y l V y 2 ) M (xyz) + tt x W (z'fi 2 - V* - z l Wl) M (xyz) ] . (43) 

Next, we commute (J42)) three times with 7r and commute (|43|) once with n and twice with 
ip to transform these operator identities into coupled differential equations for M and N: 

[x> 2 - V x 2 ) - V x ] iV x(yz) + [y\^ - V 2 ) - V y ] iV y(xz) + [z^ 2 - V 2 ) - Vj] iV z(xy) 

= -6zM(x-y)<J(x-z), (44) 

^iV y(xz) + yW. (xy) = 3 [x\^ 2 - V 2 ) - V x ] M (xyz) . (45) 

Equation is similar to Eq. (71) in Ref. Q: 

(/, 2 - V 2 )iV x(yz) + (/i 2 - V 2 )iV y(xz) + - V 2 )iV z(xy) = -65(x - y)(5(x - z). (46) 

By permuting x, y and z, we rewrite Eq. (72) of Ref. p] as 

iV y (xz) + iV z(xy) = 3(/i 2 - V 2 )M (xyz) . (47) 

This equation is similar to Eq. (|45j). The solutions for M and iV are given in Eqs. (83) and 
(84) of Ref. [J: 

4 /Y e *( x -y)-p+*( x - z )-q 
M (xy Z ) = -(2^« i/ rfprfq P (p,q) ' (48) 



where P(p, q) = 4[p 2 q 2 - (p ■ q) 2 ] + 4/i 2 (p 2 + p • q + q 2 ) + 3/i 4 , and 
iV x(yz) = 3 (V y • V z + i/x 2 ) M (xyz) 



12 /Y (-P -Q+ lu 2 ) e i ( x - y ) p+i ( x - z > q , 
Mrfprfq v — ^-f r . (49) 



{2-nfJJ *- " D(p,q) 

If these expressions for M and AT solve (|44|) and ()45|) . then we can claim that the operator 
C transforms as a scalar. To show that this is indeed true, we multiply (j4T)j) by x % and then 
subtract the result from (14411 to obtain 



(^-x i )(// 2 -v y 2 )iv y(xz) + (^-^)( / i 2 -v z 2 )iv z(xy) = v^iv x(yz) + v;iv y(xz) + v^ z(xy) . (50) 

By making the change of variables p — > q, q — > — p — q, we obtain A" y ( xz ) in the form 

12 ff, , [q-(p + q) + ±/< 2 ]^ y>P+ ' (l ~ a, ' q 

N *™ = -w i JJ dp * 1 ' ( ' 

where we have used £>(q, — p — q) = T>(p, q). Similarly, by making the change of variables 
p — > — p — q, q — > p and using the identity T>(— p — q, p) = P(p, q), we obtain 



12 /•/■ [p ■ (p + q) + \A e ^-y)-P+Kx-z). q 



Thus, from (|4*9^). (joTj) . and (|5^j) we find that the right-hand side of (|5Uj) becomes 

RW o 12 * /"/", , [p^q 2 + 2p ■ q) + q*(p 2 + 2p • q)] g^^— W 
RHS ° f (50)= (2^i/ rfprfq P(p^) • (53) 

Also, substituting (01) and 452) into the LHS of (JSUj), we get 

LHS of (50) = - JjL jfjfdpdq {( y < - x i) [ q • (p + q) + I^ 2 ] (p 2 + 

i(x-y)-p+i(x-z)-q 

- X*) [p • (p + q) + |/i 2 ] (q 2 + /i 2 )} . (54) 

We must now show that the RHS of (foTIj) in (|53~|) and the LHS of ([5U|) in (jo^|) are equal. 
To do so, we substitute the identities {y i - x i) e i^-y)-p+i(^)-i = ^ V ^e i(x - y) p+i(x - z) q and 
(z { - ^) e *(x-y)-P+i(x-z) q = iV ^( x -y). P +i( x - z ).q into ^ an d obtain 

LHS of (50) = // rf p rf q {[ q - (P + t:^ (p2 + ^V; 



(2vr)6 77 ^ ^ ^ P(p,q) 

[p-(p + q) + |/i 2 ] (q 2 + /* 2 ) ^ 1 ,( X _ y ).p +4(x - Z ). q r ,,, 

P(p,q) Vq • (55) 
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Integration by parts then yields 



LHS of (50) ^/A PC / q iv^ + q) + ^ ](p2 + ^ 



(56) 



Jp-(p + q) + ^ 2 ] (q 2 + /^ 2 ) l pl( x- y ).p +i (x^). q 

+ q 2?(p,q) J 6 

Finally, we substitute the algebraic identity 

t [q ■ (p + q) + |/i 2 ] (p 2 + /i 2 ) Jp • (p + q) + |/i 2 ] (q 2 + /i 2 ) 
X?(p,q) X?(p,q) 
= [pHq 2 + 2p-q) + q i (p 2 + 2p-q)] 

£>(p,q) 

to establish that the LHS and the RHS of (|5U|) are equal. By combining this result with 
(JUJ), we show that flSJ) holds. 

Using the same technique, we can prove that the Eq. holds as well. First, we multiply 
(JH|) by and subtract the result from (pf5j) to get 



^ - x*)iV y(xz) + (</ _ x*)iV z{xy) = -3V x M (xyz) . 



(57) 



Second, following the technique we used to derive ([51)]). we find that the LHS of (|57ji becomes 



LHS of (57) 



12i 



c/p <iq 



.i q • (p + q) + \v 2 



1 ,,2 



p • (p + q) + o/ 2 



Z>(p,q) 

i(x-y)-p+i(x-z)-q 



P 2>(p,q) 

and by substituting into (f57|). we find that the RHS of (J57)l becomes 



(58) 



RHS of (57) 



12i /"/" (p l + q i W( x ~ y )' p+ *( x ~ z )' q 



<ip dq ■ 



{2-nYJJ f n X?(p,q) 
Finally, we substitute the algebraic identity 

'" 2 p-(p + q) + V p' + q' 



(59) 



q- (p + q) + 2^_ + V ir ^ -i- 2' 



q D(p,q) p X>(p,q) P(p,q)' 

to establish that the LHS of (|57j) and the RHS of (|57j) are equal. Thus, (jBTj) holds and we 
may combine this result with P7)l to establish that (jjE]) is valid. We conclude that to order 
0(e 2 ) the operator C transforms as a Lorentz scalar. 
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IV. FINAL REMARKS 



We have shown that the operator C in a toy model iip quantum field theory transforms 
as a Lorentz scalar and that to order 0(e 2 ) the C operator in an ie(p 3 quantum field theory 
transforms as a Lorentz scalar. Based on this work, we conjecture that in general the 
C operator for any unbroken PT-symmetric quantum field theory always transforms as a 
Lorentz scalar. 

By establishing the Lorentz transformation properties of C, we can now begin to un- 
derstand the role played by this mysterious operator: Apparently, the C operator is the 
non-Hermitian PT-symmetric analog of the intrinsic parity operator Vj. A conventional 
Hermitian quantum field theory does not have a C operator. However, if we begin with 
the Hermitian theory corresponding to e = and turn on e, then the (scalar) Vj sym- 
metry of Hermitian theory disappears and is replaced by the (scalar) C symmetry of the 
non-Hermitian theory. 

The complex cubic quantum field theory discussed here is especially important because 
it contra 
gularity 

quantum field theory on a Hilbert space having a Lorentz invariant inner product. 



s the dynamics of Reggeon field theory ^| and describes the Lee- Yang edge sin- 
I5I ]. The work in this paper shows that an iip 3 field theory is a consistent unitary 
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